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Abstract. Observations at mid-infrared (4.8-17.65 /im) and radio (0.7-1.3 cm) wavelengths are used to constrain 
the structure of the high-mass star- forming region W3 IRS5 on (f.'l (200 AU) scales. Two bright mid-infrared 
sources are detected, as well as diffuse emission. The bright sources have associated compact radio emission and 
probably are young high-mass stars. The measured sizes and estimated temperatures indicate that these sources 
together can supply the observed far-infrared luminosity. However, an optically thick radio source with a possible 
mid-infrared counterpart may also contribute significant luminosity; if so, it must be extremely deeply embedded. 
The infrared colour temperatures of 350-390 K and low radio brightness suggest gravitational confinement of the 
H II regions and ongoing accretion at a rate of a few 10~ 8 Mq yr _1 or more. Variations in the accretion rate 
would explain the observed radio variability. The low estimated foreground extinction suggests the existence of a 
cavity around the central stars, perhaps blown by stellar winds. At least three radio sources without mid-infrared 
counterparts appear to show proper motions of ~100 km s _1 , and may be deeply embedded young runaway OB 
stars, but more likely are clumps in the ambient material which are shock-ionized by the OB star winds. 

Key words. Stars: Circumstellar matter; Stars: formation; Instrumentation: high angular resolution 



1. Introduction 

Stars of masses >8 M Q spend a significant fraction of 
their lifetimes, >10%, embedded in their natal molecular 
clouds. Single-dish (sub)millimeter observations have clar- 
ified the structure of hi gh-mass protostellar envelopes on 
-10 4 -10 5 AU scales (see lHatchell k van der Takl2003l and 
references therein). However, the distribution and kine- 
matics of material on < 1000 AU scales is poorly known, 
due to the large ( > 1 kpc) distances involved, and the 
lack of tracers at optical and near-infrared wavelengths. 
These scales are of great interest to decide between for- 
mation mechanisms for high-mass stars, an d to clarify 
the r elation with clustered star formation l|ChurchweJ 
120021) . Does the distribution of stellar masses in a star- 
forming region depend on the stellar density? Also, the 
origin of the observed outflows and their interaction 
with the environment on < 1000 AU scales remain un- 
clear. Subarcsecond resolution observations are necessary 
to shed light on these and other questions, which, at 
(sub )millimcte r wave lengths, are just coming within reach 
(Beu ther et al.ll2004[) . However, these resolutions can al- 



ready be achieved in both the infrared and radio wave- 
bands, where extinction is much smaller than in the opti- 
cal. 

In the infrared, high-resolution techniques are most 
advanced at near- infrared wavelengths. Such observations 
probe less embedded, more evolved phases, where a signifi- 
cant part of the surroundings is already ionized. Important 
progress has been made with the identification of the 
ionizing stars of several ultracompact H II regions (e.g., 
IWatson k Hansonl Il997l iFeldt et aH 12003^ . In the case 
of intermediate-mass stars, the imaging of the hot inner 
regions of disks is presently generatin g a lot of interest 
l|Danchi et all200 ltlTuthill et al.l20fJl|) . In addition, a few 
more embeddded objects ha ve been probed l|Weigelt et alJ 
l2002t IPreibisch et aljEtoi . although in the general case 
long baseline interferometers will be needed to tackle 
most targets given the characte ristic size scales involved 
llMonnier k Millan-Gabet] l2002|) . At mid-in frared wave- 
length s, pioneering work has been done by I Walsh et alJ 
(2001), but s ubarcsecond resolution has only recently 
been achieved jTuthill et alJl2002t iDe Buizer et alJl2002t 
iGreenhill et al.ll2004 - 
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In the cm-wave region, most subarcsecond-resolution 
studies have concentrated on H2O masers, which are 
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Fig. 1. Images of W3 IRS5 at 43 GHz (left) and 22 GHz (right). Contours start at 4cr and increase by 2cr, where 
er=0.16 mJy/beam at 43 GHz and 0.15 mJy/beam at 22 GHz. 



bright enough for Very Long Baseline Interferometry 
(VLBI) observations. The excitation requirements of the 
masers are such that the emission usually traces shocks as- 
sociated with infalling or outflowing motions. The VLBI 
data indicate that the maser emission traces moving gas 
parcels, rather than shock waves propagating through an 
H2O rich cloud. In the case of outflow motions, both bipo- 
lar and spherical flows are seen, which may rep resent dif- 
ferent stages of evolution ijTorrelles et al.ll2003[) . In some 
cases, H2O masers in star-formi ng; regions may arise in 
accre tion shocks in infalling gas l|Menten & van der Takl 
120041) . 

Continuum emission at centimeter wavelengths arises 
in ionized gas. In stellar winds and outflows, the gas can 
be collisionally ionized, and VLBI data indicate a mixture 
of bipolar and equatorial outflows llHoarell2002<) . Close to 
hot stars, small regions of photo-ionized gas are observed 
as 'hypercompact' H II regions, which represent a very 
early stage of high-mass star formation ( Garav fc Lizanol 

[3. 

1.83±0.14 kpc (|lmai et alJ 
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Orion. The brigh t mid-infrared sourc e has been resolved 
into a double l|Howell et al.l Il98li iNeugebauer et all 
Il982h . Single-dish submillimeter mapping indicate an 
envelope mass of 262 M Q within a ra dius of 60,000 AU, 
with an r -15 density distribution l|van der Tak et all 
12000^1 . Near-i nfrared imaging sh ows a dense cluster 
(~3000 pc~ 3 : iMegeath et al.lfl996T) . mostly composed of 
low-mass pre- ma in- sequence stars with ages 0.3-1 Myr 
ijOiha et al.l E004L Radio continuum observations show a 
cluster of at least six 'hyperc o mpact' H II regions lab eled 
A...F l|Claussen et alJ Il994t iTieftrunk et"afl Il99it . at 



least one of which exhibits proper motions l( Wilson et alJ 
120031) . Water maser mapping reveals ^100 spots, grouped 
in two flows: one roughly spherical and centered close 
to continuum source A, and the other more collimated 
and centered cl ose to source D ( Claussen et al.1 119941 
llmai et al . 2000). Mid-infrared spectroscopy shows CO 
absorption features blueshi fted by 4.. 46 km s -1 relative 
to the systemic velocity ijMitchell et alJ Il99l|) . which 
must arise in an outflow. In millimeter-wave CO emission, 
blue- and redshifted outflow lobes are detected out to 
«23 km s _1 from the systemic velocity I Claussen et alJ 
Il984h . indicating that the highest-velocity gas is very 
compa ct. Finally, Chandra observations bv iHofner et alJ 
(2002) indicate an X-ray luminosity of W3 IRS5 of 
£x=9xl0 29 erg s _1 (for d=1.83 kpc), consistent with the 
typical values for T Tauri stars. 

This paper presents new cm- wave and mid-infrared im- 
ages of W3 IRS5 at sub-arcsecond resolution. The goals 
are to clarify the nature of the radio continuum sources 
and their relation with the infrared double, to find which 
ones are self-luminous, and which ones power the region. 

2. Observations 

2.1. Radio Observations 

Radio observations of W3 IRS5 were carried out with 
the NRAO 1 Very Large Array (VLA) on 1996 October 
24, when the VLA was in its A-configuration. At this 
time, thirteen VLA antennas were equipped with 43 GHz 
receivers; the other fourteen observed at 22 GHz (re- 
spectively known as Q- and if-band in radio astron- 



1 The National Radio Astronomy Observatory (NRAO) is 
operated by Associated Universities, Inc., under a cooperative 
agreement with the National Science Foundation. 
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Fig. 2. Mid-infrared long-exposure images of W3 IRS5. Contours are at 0.5, 1, 2, 3, 5, 10, 20, 30 and 70% of the peak 
intensity. 



omy). Zenith opacity was 0.089 at 22.5 GHz and 0.071 
at 43.3 GHz. Elevation-dependent antenna gains were in- 
terpolated from values measured by the VLA staff. The 
phase calibrator, 0228+673, was observed every 10 min- 
utes at 43 GHz and every 15 minutes at 22 GHz (the 'fast 
switching' procedure was not implemented at the time). 
Pointing was checked at 8.4 GHz on the same source ev- 
ery 70 minutes at 43 GHz and every 5 hours at 22 GHz. 
On-source integration time is 438 min at 22 GHz and 
400 min at 43 GHz. The data were edited, calibrated, 
and imaged with NRAO's Astronomical Image Processing 
System (AIPS). Absolute calibration was obtained from 
observations of 3C 286 using flux densit ies interpolated 
from the values given bv lOtt et ail l|l994|) . For 0228+673, 
we obtain a flux density of 1.83 Jy at 22 GHz and 1.55 Jy 
at 43 GHz. 

Figure H shows the 22 and 43 GHz maps, which have 
rms noise levels of 0.15 and 0.16 mJy beam -1 . These 



maps were obtained from the uv data by a Fourier trans- 
form with uniform weighting, and deconvolved with the 
Clean algorithm. Restoring beam major and minor axes 
and position angles are 44x37 mas at position angle -68° 
at 43 GHz and 89x88 mas, PA -53° at 22 GHz. 

2.2. Infrared Observations 

Data were obtained in August 2002 with the Long 
Wavelength Spectrometer (LWS) camera on the Keck I 
telescope 2 . Two different observing methodologies were 
employed. The first set of observations utilized a standard 



2 The W. M. Keck Observatory was made possible by the 
support of the W. M. Keck Foundation, and is operated 
as a scientific partnership among the California Institute of 
Technology, the University of California, and the National 
Aeronautics and Space Administration. 
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Fig. 3. Mid-infrared 'speckle' images of W3 IRS5. Contours are at 0.5, 1, 2, 3, 5, 10, 20, 30 and 70% of the peak 
intensity. 



chop-nod pattern, with frames coadded to build up longer 
exposure times. Observations of W3 IRS5 (with the cali- 
brator £ Cyg) taken in this mode emphasized the recovery 
of faint structure. The second mode had a fast readout in 
addition to the chop-nod, so that large volumes of rapid 
exposure data were collected in a fashion analogous to a 
mid-IR 'speckle' experiment, with the hope to recover fine 
structure in the images. Again, these objects were paired 
with similar data taken on a point source calibrator, this 
time a Ceti. No point-source calibrator files were taken 
for 9.9 /xm or 12.5 /im observations. These two modes are 
denoted 'L' and 'S' for 'long' and 'speckle' exposure here- 
after. 

Data have been analyzed with an iterative matched fil- 
ter version of the shift-and-add algorithm. This algorithm 
attempts to match the shifts in the current iteration to 
maximise the correlation with the output of the previ- 
ous iteration. Significant gains in image resolution were 



demonstrated over straight coadding, while a simple shift- 
and-add strategy was foiled in this case by the presence 
of two nearly equal peaks. 

Figures |21 and |3| present the resultant images, except 
the 17.65 /Ltm long-exposure image which is very similar 
to the speckle image at that wavelength. Note that an ar- 
tifact due to an unwanted reflection from an optical sur- 
face affects the 4.7 /jm image, giving a spurious feature 
to the North-East which is not seen at any other wave- 
length. This 'ghost' was also present in point-source cali- 
brator data in this filter. The 9.9 [im and (to a lesser ex- 
tent) the 10.7 /jm data have degraded signal-to-noise due 
to the much lower flux levels. There were, unfortunately, 
additional experimental difficulties which were not easy 
to account for. The observations were taken under condi- 
tions of variable cirrus, increasing the errors in photom- 
etry. Compounding this was an intermittent mechanical 
fault with the camera mechanism which resulted in par- 
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tial occultation of the pupil, affecting both the throughput 
and the point-spread function (PSF). Although this had 
little effect on the maps presented here, it did preclude 
our original intent of fully deconvolving the images using 
the PSF from the reference star observations. 



3. Results 

3.1. Radio Positions 

Table □ lists the sources detected above 5<r using the 
multiple-peak-finding A1PS task SAD after trying various 
weighting schemes. For the Q-band data, the table uses 
the image obtained with uniform weighting, which has bet- 
ter positional accuracy while for the i^-band data, tabu- 
lated results use an image obtained with natural weight- 
ing, where more sources are detected (rms = 142 /iJy). 
Sources Ql and Q4 are 4a detections and were not found 
by SAD, but their detection is secure because they have 
-fT-band counterparts. The flux densities in the table are 
corrected for primary beam response. Positional uncer- 
tainties are statistical errors and apply to the relative po- 
sitions of these radio sources. 

Columns 6-8 of Table gives the source sizes. Since 
only sources Q3, K6 and K8 appear marginally resolved, 
the sizes have not been deconvolved. The sizes would be 
lower limits i f extended emission is resolved out by the in- 
terferometer l(Kurtz et al.lll 999: Ki m fc Kooll200ll) . In the 
particular case of W 3 IRS5, however, multi-configuration 
observ ations rule out extend ed emission down to very low 
limits (Tieftrunk et al .119971) . Therefore, the radio sources 
of W3 I RS5 belong to the cla ss of 'hypercompact' H II re- 
gions l|Kurtz fe Francoll2002|) . 

C omparing the positions in Table ^ wi th th ose from 
1989 l|Claussen et"alll994 iTieftrunk et aTlll997|) leads to 
the following identifications: Ql = K2, Q2 = K3 = A, Q3 
= K4 = B, Q4 = K6, Q5 = K7 = MD1, K5 = C, K8 
= F. Source Kl is several arc seconds away and proba- 
bly unrelated. Sources Q1=K2 and Q4=K6 were not seen 
before and appear to be new. On the other hand, sources 
E and G seem to have disappeared since 1989, the epoch 
when the Tieftrunk et al data were taken. Most strikingly, 
source D2 has disappeared, which is remarkable since it 
was the strongest source in 1989. Perhaps Dl and D2 were 
not separate sources, but merely substructure within one 
source, which we refer to as D hereafter. 

Our sources with counterparts in the old data do not 
exactl y lie on the position s repo rted by Idaussen et alJ 
()l994f) and iTieftrunk et all l)l997j) . but rather at 80-140 
mas shifts. The shifts are 2-3 beam sizes and >10 times 
the formal error on relative positions. The position an- 
gles of the shifts vary between 20° and 60°, which ar- 
gues against instrumental effects such as pointing er- 
rors or changes in calibrator positions, which would shift 
all sources in the same direction. The data thus seem 
to confirm the exis tence of proper motions reported by 
IWilson et al.l l)2003h . At a distance of 1.83 kpc, a motion 
of 100 mas in 7.48 yr corresponds to a transverse velocity 



of 116 km s _1 . The space velocity may be a factor y/2 
higher, or 164 km s _1 . These values are much larger than 
the m otions of the H2O masers of ~20 km s _1 ijlmai et alJ 
2000) and of the CO emission and absorption (§0. 



3.2. Infrared Positions 

The infrared images (Figs. [5] and 03 show three compact 
sources, surrounded by diffuse emission which becomes 
more pronounced toward longer wavelengths. We begin 
a quantitative analysis of these images by fitting simple 
profiles to the data, and by measuring flux densities in 
different regions. The results are given in Table which 
gives the fluxes, relative positions and full-width at half- 
maximum (FWHM) of the various components. In addi- 
tion, the FWHM of the point-source reference stars are 
given, which gives an estimate of the resultant system 
PSF. Examination of these data shows that the system 
appears truly diffraction-limited in either mode ('L' or 'S') 
at 17.65 /mi. At shorter wavelengths (8.0-12.5 /um), the 
'S' mode delivers a significantly smaller FWHM than 'L' 
(8.0 /jm gives a direct comparison) which approaches the 
formal diffraction limit. The dramatic increase in size at 
4.8 [im implies some optical problem beyond the normal 
effects of diffraction and seeing, such as optical aberration 
or defocus. 

In the absence of a wide-field image with standard 
stars, our only astrometric information comes from the rel- 
ative positions of the components. We refer to the north- 
ernmost bright component as MIR1, with MIR2 being of 
nearly equal brightness to the south, and MIR3 the much 
fainter southernmost peak. Table lists the separation 
and the position angle of MIR2 & MIR3 relative to MIR1 
for all observations, obtained through Gaussian fits to the 
emission. 

The mean separation of MIR1 & MIR2 is 1124±74 
mas at a position angle of 36.8±1.7 degrees. These values 
are consistent with those from earlier mid -infrared work 
jHowell et al.lll98lilNeugebauer et al.ll982j) . but our data 
are the first to image the mid-infrared double directly. 

We compare this relative position with those of pairs 
of radio sources. The best match is for pair Q3-Q5, whose 
separation of 1210 mas at a position angle of 37?4 is in 
good agreement with the infrared peaks. The only other 
radio pair with similar relative positions is K7-K8, which 
has a separation of 1141 mas, but at a position angle of 
28?9, inconsistent with the infrared result. On this ba- 
sis, we identify the bright mid-infrared sources with radio 
sources Q5=K7=MIR1 and Q3=K4=MIR2. The position 
of source MIR3, with 1-2% of the flux density of the main 
sources, then coincides with radio source Q4=K6, con- 
firming our identification. There is a cluster of H2O maser 
spots close to this obj ect, at A(a)«250 and A((5)w2000 
mas l|Imai et a 
identifications. 



ti is obi e 
. 2000). 



Table summarizes our source 



The images in Figures and |21 show the only region of 
flux detected within the 10'.'24 field of view of LWS, with 
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Table 1. Radio emission from W3 IRS5. Numbers in brackets are uncertainties in units of the last decimal. 



Source 


a 


5 


Peak I v 


Total S v 


Major axis 


Minor axis 


PA 




(J2000) 


(J2000) 


mjy/beam 


mjy 


mas 


mas 


deg 


Ql 


02 25 40.660407(635) 


62 05 51.82215(404) 


0.691(162) 


0.78(31) 


45(10) 


41(10) 


93(90) 


Q2 


02 25 40.676344(379) 


62 05 52.04937(261) 


1.285(158) 


1.92(36) 


59(7) 


41(5) 


47(13) 


Q3 


02 25 40.681511(211) 


62 05 51.53000(198) 


2.102(155) 


3.89(42) 


63(5) 


47(3) 


7(10) 


Q4 


02 25 40.728475(704) 


62 05 49.85180(702) 


0.666(154) 


1.43(46) 


75(17) 


46(11) 


21(18) 


Q5 


02 25 40.783441(212) 


62 05 52.46552(141) 


1.996(163) 


2.19(30) 


43(4) 


41(3) 


88(57) 


Kl 


02 25 40.143241(1234) 


62 05 52.28144(1080) 


0.745(142) 


0.85(27) 


139(27) 


99(19) 


25(22) 


K2 


02 25 40.661233(1029) 


62 05 51.86622(1174) 


0.847(140) 


1.08(29) 


175(29) 


88(15) 


160(9) 


K3 


02 25 40.675758(1135) 


62 05 52.04087(690) 


0.857(142) 


0.76(23) 


121(20) 


88(15) 


59(21) 


K4 


02 25 40.682172(526) 


62 05 51.54376(441) 


1.938(139) 


2.76(31) 


152(11) 


113(8) 


154(10) 


K5 


02 25 40.712459(785) 


62 05 52.44786(569) 


1.163(142) 


1.12(24) 


110(13) 


106(13) 


168(138) 


K6 


02 25 40.733047(1601) 


62 05 49.84977(1214) 


0.797(134) 


1.76(41) 


174(29) 


154(26) 


25(56) 


K7 


02 25 40.782747(737) 


62 05 52.46103(637) 


1.313(140) 


1.63(28) 


150(16) 


100(11) 


29(10) 


K8 


02 25 40.861211(888) 


62 05 53.46070(717) 


1.307(135) 


2.44(37) 


169(17) 


135(14) 


25(19) 



Table 2. Quantitative analysis of the infrared imaging: flux densities and relative positions. The table columns are 
arranged as follows: (1) observing A/AA; (2) speckle 'S' or long 'L' exposure; (3) FWHM of point-source reference 
star; (4,6,10) flux of MIR 1-3; (5,7,11) FWHM of MIR 1-3; (8,12) separation of sources MIR1-2 and MIR1-3; (9,13) 
position angle of sources MIR1-2 and MIR1-3; (14) total flux detected. 



Filter 


Mode 


PSF star 




MIR1 




MIR2 


MIR1-MIR2 


MIR3 


MIR1-MIR3 


Total 


A/AA 


L/S 


FWHM 


flux 


FWHM 


flux 


FWHM 


Sep 


PA 


flux 


FWHM 


Sep 


PA 


Flux 


(H 




(mas) 


(Jy) 


(mas) 


(Jy) 


(mas) 


(mas) 


(deg) 


(Jy) 


(mas) 


(mas) 


(deg) 


(Jy) 


4.8/0.6 


L 


507 


44 


453 


23 


471 


1217 


215 


0.6 


480 


2730 


186 


80 


8.0/0.7 


L 


384 


140 


432 


108 


480 


1198 


215 


2.0 


436 


2741 


187 


301 


8.0/0.7 


S 


258 


152 


310 


123 


367 


1214 


217 


1.8 


304 


2663 


188 


328 


9.9/0.8 


L 




4 


576 


4 


568 


1013 


215 


< 0.5 








11 


10.7/1.4 


S 


291 


13 


460 


12 


552 


1078 


219 


< 1.2 








32 


11.7/1.0 


L 


364 


60 


519 


50 


607 


1097 


216 


0.7 


255 


2739 


187 


142 


12.5/0.9 


S 




60 


467 


50 


587 


1145 


218 


1.1 


329 


2632 


188 


146 


17.65/0.9 


L 


451 


98 


716 


92 


636 


1049 


217 


< 10.7 








318 


17.65/0.9 


S 


452 


107 


646 


93 


615 


1102 


219 


< 10.3 








323 



one exception. In a few frames (which happened to be off- 
set from center) a faint diffuse component was seen at the 
extreme edge of the field, 7'.'3 from MIR1 at a position 
angle of 160°. Using the radio identifications of MIR1 and 
MIR2 as astrometric reference, this source, which we call 
MIR4, lies at position a = 02' 1 25 m 41!1388, 6 = 62° 05' 
45'.'604 (J2000), where no radio emission is detected. The 
extended nature and location at the edge of the field of 
view preclude measurement of its mid-infrared flux den- 
sity. 

3.3. Radio Brightness 

The flux densities of 22 GHz sources K3, K4, K5, 
K7 and K8 are si g nifica ntly different from the values 
by iTieftrunk et al. I <1997t> . probably due to variability. 
Instrumental effects, such as calibration problems, atmo- 
spheric decorrelation, or difference in beam size, would 
affect all sources in the same way. Instead, several com- 
ponents seem to undergo gradual increases or decreases 
in 22 GHz brightness over the available 13-year period 
(Figure EJ. 



Table 3. Radio and infrared identifications. 



MIR 


Q 


K 


1989 


1 


5 


7 


D 


2 


3 


4 


B 


3 


4 


6 






1 


2 






2 


3 


A 






5 


C 






8 


F 








E 








G 






1 





The total flux densities in Table Q indicate a spectral 
index 7, defined through 

S v cx v 1 

of 7=1.42±0.77 for source Q2, 7=0.52±0.33 for source Q3, 
and 7=0.45±0.48 for source Q5. Values derived from the 
peak brightness are 0.62±0.44, 0.12±0.22 and 0.64±0.29. 
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These values are consistent with thermal emission, and 
are not affected by variability as the data were taken si- 
multaneously. For the weak 43 GHz sources Ql and Q4, 
we find 7=-0.3±0.6, i.e., flat or slightly nonthermal spec- 
tra. Sources Kl, K5 and K8, which are not detected at 
43 GHz, have 7=-l±l and may be of nonthermal nature. 
T he spectral ind i ces fo und here are consistent with those 
bv lWilson et alJ l|2003|) . which are also based on simulta- 
neous measurements, except for source K8=F which was 
detected at 43 GHz in 2002 but not in 1996. 

In the case of Bremsstrahlung from an ionized region 
with a power law distribution of the electron density with 
radius, 

-<? 



n e oc r 



the spectral index is 



7 =(2< Z -3.1)/( (? -0.5) 



ll01nonlll97,^ . Hence the spectral index of 7 « 0.5 mea- 
sured for Q3 and Q5 corresponds to q 1.9, close 
to the value of 2 for an ionized wind. Observations of 
broad H I radio rec ombination lines ( Tieftrunk et al.ll997t 
ISewilo et al1l2004h support this interpretation, although 
higher angular resolution is needed to relate the line emis- 
sion with the continuum sources. In any case, our mea- 
surements are not sensitive enough to rule out a flat radio 
spectrum for these sources, which would indicate optically 
thin emission. The objects could then be externally ion- 
ized, consistent with the absence of mid-infrared emission. 

The value 7 w 1.5 measured for Q2 could arise in an 
H II region with a constant density at the center and a 
steep (q sa 4.7) outer fall-off. It is not clear which mecha- 
nism would create such a density distribution, especially 
since the sound speed of vs ~10 km s _1 of H II regions im- 
plies that in these compact (0 < 100 AU) sources, density 
fluctuations are washed out within ~50 yr. Therefore, Q2 
is probably a uniform-density H II region which is (moder- 
ately) optically thick, but again, the data do not rule out 
a wind spectrum. In either case, it is internally ionized. 

More sensitive measurements over a larger wavelength 
range are necessary to constrain the emission mechanism. 
A constant spectral index would support the wind model, 
while a bent spectrum would indicate uniform H II regions 
of intermediate optical depth. Care has to be taken, how- 
ever, not to include dust emission at hig h frequencies, or 
synchrotron emi ssion at low frequencies IjFelli et alJ ll993: 
iReid et al.lli"995i 

3.4. Infrared Brightness 

Columns 4, 6 and 10 of TableElreport the flux densities of 
the three mid- infrared sources. Brightness was measured 
in circles of radius 600 mas, which cover all of the diffrac- 
tion and seeing patterns. For the purposes of flux cali- 
bration, six reference stars were used, and four of these 
gave consistent photometric results (the other two, pre- 
sumably affected by non-photometric conditions, were ig- 
nored). At 9.9 and 12.5 /mi where there were no calibrator 



star data taken in an identical way, the flux calibration 
was derived from indirect measurements and should be 
regarded as tentative. Due to these difficulties and the 
variable vignetting in the camera, we quote errors on the 
photometry of up to 50%. Within these errors, the total 
flux densities at 4.8-10.7 /mi ar e consi stent with the val- 
ues measured by IWillner et al.l l|l982() in a ~10" beam. 
However, at 9.9 and 12.5 /im, where no calibrators were 
observed, and at 11.7 /mi, where only one calibrator was 
observed, the photometric error is probably closer to a 
factor of 2. Indeed, at 11.7 and 12.5 /im, our flux densities 
are a factor of ~2 below the v alues measured by Willncr 
et al and bv lPersi et al] l|l996|) in a 3" beam. 

It is interesting to note that the relative fluxes and 
positional offsets between MIR1-3 remain fairly constant 
across the mid-infrared and (presuming our radio identi- 
fications) into the radio. This implies that it is unlikely 
that there are large differences in the effective tempera- 
tures or the optical depths to these 3 components. The 
only readily identifiable exception to this is a systematic 
trend for MIR1 being brighter than MIR2 at short wave- 
lengths, while they are nearly equal at long wavelengths. 
This may imply a somewhat hotter underlying spectrum, 
or there may be opacity gradients in the line of sight. 

3.5. Contribution from PAHs 

The observed infrared emission may be continuum emis- 
sion from dust grains. However, the mid-infrared spectra 
of many Galactic sources, including compact H II regions 
and other star-forming regions, show strong emission fea- 
tures due to Polycyclic Aromatic Hydrocarbons (PAHs). 
The stronges t PAH features lie a t 3.3, 6.2, 7.7, 8.6, 11.2 
and 12.7 /mi 1 Peeters et al.ll2002j) . Therefore our 8 /mi fil- 
ter contains the 7.7 /mi feature, the 10.7 and 11.7 /mi 
filters the 11.2 /mi feature, and the 12.5 /mi filter the 
12.7 /mi feature. The strength of these features reflect the 
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Fig. 4. Flux densities (at 22 GHz) of radi o sources in W3 
IRS 5 as a function of time. Data are from lTieftrunk et alJ 
Jl997lUWilson et all ll2003h . and Tabled 
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ambient UV radiation field, rather than the dust mass or 
temperature. Therefore, to interpret our Keck data prop- 
erly, quantifying the contribution of PAHs to the observed 
emission is essential. 

We have searched the ISO-SWS spectrum of W3 IRS5 
(F. Lahuis, priv. comm.) for PAH features. With typi- 
cal widths of 0.1-0.4 fxm, the features should be easily 
resolved with ISO. No PAH features are detected down 
to an rms noise level of Jy. The flux densities mea- 
sured with ISO are 60-90% of those measured with Keck, 
so beam dilution does not play a role. We conclude that 
the emission observed with Keck is thermal emission from 
dust grains. 

3.6. Infrared Sizes 

Our fits to the mid-infrared images show that the two 
bright cores MIR1 and MIR2 exhibit systematically larger 
sizes than the reference stars, as measured by the Gaussian 
FWHM values given in Column 3 of Table |21 In this sec- 
tion, we extract quantitative estimates of the apparent an- 
gular diameters of these cores, by deconvolving with the 
reference star PSF then fitting with a simple circularly- 
symmetric profile (in this case a uniform disk). 

However, this could only be done in a minority of cases 
where the data were suitable and of sufficient quality. The 
4.8 /xm data were affected by an unknown optical prob- 
lem (as discussed earlier), while the 9.9 & 12.5 /im data 
had no PSF reference star data, and were ignored here. 
Furthermore, deconvolution requires the highest possible 
angular resolution data, and we therefore restrict our at- 
tention to only the rapid exposure observing mode 'S', dis- 
carding 'L' (e.g. all 11.7 /im data). Perhaps the most diffi- 
cult aspect of the deconvolution problem was distinguish- 
ing between the resolved cores and the more extended neb- 
ula. At the longest wavelength, 17.65 /im, this was not 
possible for two reasons: firstly the extended component 
was relatively bright compared to the cores, and secondly 
the angular resolution was not sufficient to clearly distin- 
guish between compact and extended flux. 

The remaining datasets suitable for deconvolution and 
diameter fitting were from 8.0 & 10.7 /im. Fits were ob- 
tained with a uniform circular disk profile, although in the 
partially resolved case as here any simple model (such as a 
Gaussian) would serve equally well. Uniform disk fits were 
obtained for MIR1 & MIR2 at 8.0 /im, where the diame- 
ters are 207 & 254 mas, and at 10.7 /mi, where they are 
300 & 333 mas, respectively. Although the formal errors 
on these quantities are around 40 mas, the true uncer- 
tainties are hard to quantify due to unknown seeing and 
optical changes between the source and calibrator stars, 
and due to imperfect rejection of the extended nebula in 
fitting the core. 

Two systematic trends are noted here: MIR2 appears 
slightly larger than MIR1, and the sizes at 10.7 /jm are 
larger than those at 8.0 /im. However, particular caution 
needs to be expressed over contamination from the ex- 
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Fig. 5. Mid-infrared spectral energy distributions of W3 
IRS5 after dereddening. Squares, triangles and dots in- 
dicate MIR1, MIR2 and MIR3. Model curves are super- 
posed. 

tended flux component which may have a role in causing 
these apparent extensions. 

The measured mid-infrared sizes exceed the limits on 
the radio sizes of ~100 AU (Tabled]). This result supports 
a model where the radio emission comes from ionized gas 
very close to a star and the mid-infrared emission from 
warm dust somewhat further out. 

4. Discussion 

Our observational findings can be summarized as follows: 
three mid- infrared point sources with radio counterparts; 
four radio sources without mid-infrared counterparts that 
appear to change position; and diffuse mid-infrared emis- 
sion. The following sections discuss each of these compo- 
nents in turn. 

4.1. Mid-infrared point sources 

To estimate the physical properties of the compact mid- 
infrared sources, we have compared their flux densities to a 
simple blackbody model. Given the lack of strong observed 
colour variations (§ I3.4fl . we assume that the sources have 
the same temperature and foreground extinction. Based 
on the results of § 13.61 we use a radius of 250 AU for sources 
MIR1 and MIR2; for MIR3, i?=30 AU is adopted based on 
its lower brightness. The observed far-infrared luminosity 
then limits the temperature to T <390 K, which appears 
plausible base d on the measured colours and the 2.2 /tm 
photometry bv lOiha et ail l|2004[h 

Using this temperature and radius, we can model the 
observed flux densities if we know the foreground extinc- 
tion. The broad-ban d mid-infrared spectrum presented by 
IWillner et al.l l|l982F) indicates a silicate optical depth of 
ts «4. 3-5.0 assuming pure absorption, or ts=7.64 when 
correcting for underlying silicate emission, as Willner et 
al do. More recent data from ISO give consistent results 
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(F. Lahuis, priv. comm.), even though they refer to a 
larger beam (~20") that varies by a factor ^4 over this 
wavelength range. 

Figure [S] shows the measured flux density spectrum af- 
ter de-reddening by rs=5.0. Values for the extinction at 
other wavelengths are compu ted using dust properties by 
lOssenkopf fc Hennind l|l994|) . Model 5. This dust model 
gives a good match between envelope m asses derived from 
dust and CO ijvan der Tak et al.lll999h . The figure shows 
that a blackbody model with T=390 K reproduces the 
data within a factor of 3. The largest outliers are the 
9.9 and 11.7 /ini points, at which wavelengths the cali- 
bration is the most uncertain. 

The lOssenkopf fe Hennind l|l994j) dust mo del has 
Si:C= 1.45, while values up to «2 are observed ^Krugel 
Increasing the Si:C ratio would improve the match 
between data and model at 9.9 /im, but would give a worse 
fit at 11.7 /im. More likely, the assumption of blackbody 
emission is not quite valid, so that towards shorter wave- 
lengths, smaller radii and higher temperatures are probed. 
Geometrical effects may also influence the shape of the sil- 
icate absorption. 

Temperatures below 390 K are energetically allowed, 
but require lower extinctions to fit the observed flux den- 
sities. For T < 350 K, the required extinction drops 
belo w t,s i =4.3, which w e consider unlikely based on 
the IWillner et alJ l)l982l) data. The total luminosity for 
T=350 K is ~8xl0 4 L , which leaves -4xl0 4 L for 
a third power source, such as radio source Q2 (§ 13. 3|) . 
This option is more likely than the case of two power 
sources, because radio sources Q2, Q3 and Q5 are of sim- 
ilar strength (§ 13.3(1 . 

The (tentative) size increase of MIR1 & 2 from 8.0 to 
10.7 /im f§ !3.6fl is to be expected if a more realistic assump- 
tion of a centrally-heated dust cloud with a thermal profile 
is adopted, rather than a blackbody at a single temper- 
ature. A detailed understanding of this deeply embedded 
and complex region will clearly require radiative transfer 
modelling, and it is encouraging that mid-infrared imaging 
appears capable of placing meaningful constraints. 

The best-fit extinction of ts=5.0 is significantly below 
Willners estimate. Perhaps their simple formula to correct 
for silicate emission is not valid at large extinction values. 
More likely, the silicate absorber is physically decoupled 
from the underlying continuum emitter. One geometrical 
interpretation is that of two star/disk systems surrounded 
by a cavity, whose walls cause the silicate absorption. Such 
a cavity would also explain why submillimeter imaging in 
a 15" beam indicates a m uch higher extinction (Ay ~300; 
van der Tak et a"i]l2000j) than the mid-infrared data. 

4.2. Stationary radio sources 

The flux densities S v of the radio sources (Table ^) can 
be used to estimate the Lyman continuum emission Nl of 
their ionizing sources, assuming t hat the H II regions ar e 
uniform and isothermal (see, e.g.. lRohlfs fc Wilsonll200oh . 



This discussion concentrates on sources Q2, Q3 and Q5 
which have positive spectral indices 13.31) . and consid- 
ers both optically thin and optically thick emission as 
limiting cases. In the optically thin case, Nl is directly 
proportional to the flux density. In the optically thick 
case, black body emission at T=10 4 K indicates radii 
of ^20 AU, consistent with the observational upper lim- 
its (Table 0. The emission measure follows from setting 
the free-free optical depth equal to unity; radii and emis- 
sion measures together indicate electron densities of 10 6 - 
10 7 cm" 3 . Balancing photo ionization with 'case B' recom- 
bination ljOsterbrocklll989|) finally gives Nl- The results 
for both cases are Nl=1...7x 10 44 s _1 , with a weak depen- 
dence on electron temperature. 

The similar values of Nl derived for 43 GHz sources 
Q2, Q3 and Q5 indicate that they have similar luminosi- 
ties. Therefore it is hard to see how only the ionizing 
source of Q2 could be invisible in the mid-infrared, unless 
it is extremely deeply embedded. In fact, our 17.65 /im 
images may show a source about 0'.'5 North of MIR2, but 
the data do not allow to extract a flux density. 

The stellar luminosities of «40,000 L Q (§ 14. 1|) corre- 
spo nd to masses of ~20 M q and ZAMS spectral types 
08 IjMaeder fc Mevnetlll989tl . Their expected Lyman con- 
tinuu m emissions are «6xl0 48 s -1 l|Schaerer fc de Koterl 
Il997l) . which is ~10 4 times the value just derived from the 
radio continuum emission. Since dust absorption inside the 
H II region only accounts for factors of 2-3, th is discrep- 
ancy may be due to accretion of dust particles ( Walmslevl 
1995J). The observed variability (Fig. may then corre- 
spond to variations in the accretion rate. 

Accretion of dust would also explain why the hyper- 
compact H II regions stay confined to a ~20 AU radius. 
The ionization front around an O-type star on th e main 
sequence is a D-critical front fe.g.. lOsterbrockll989h which 
moves at about the sound speed of «10 km s -1 , or some- 
what less (5-7 km s -1 ) since the surroun ding H I shell 
needs to be accelerated. lAcord et all ()l998|) have seen such 
expanding motions in the ultracompact H II region G5.89. 
In the case of W3 IRS5, expansion at 5-10 km s _1 would 
lead to an increase in radius from 20 to 100 AU over the 
observed 10-year period which is not observed. 

The accretion rate needed to confine the hypercompact 
H II regions may be estimated by equating the accretion 
force (momentum transfer rate) of the dust to the ther- 
mal pressure of the H II region. Using a radius of 20 AU, 
a density of 3x 10 6 cm" 3 , and T=8000 K for the H II re- 
gion, we find M=1.5xl0~ 8 M yr _1 assuming that the 
dust is in free fall onto a 20 M star. In reality, radiation 
pressure will slow down the dust from the free-fall speed 
(42.2 km s _1 ), so that perhaps twice this M is needed. 
If the stars have winds with substantial mass loss rates 
(e.g., 10~ 6 M Q yr _1 ), even higher accretion rates may be 
needed to confine the H II region. 

Recent work bv lKetol 1 2002'!. however, shows that stel- 
lar gravity prevents the hydrodynamical expansion of H II 
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regions inside a 'gravitational radius' 
r g = GM/2c 2 s 

where G is the gravitational constant, M the stellar mass, 
and c s the sound speed of H II regions of «10 km s . For 
the bright mid-infrared sources in W3 IRS5, M «20 M Q 
(§ HUt so that r g «90 AU, consistent with the observa- 
tional limits (Table [TJ. 

In Keto's model, both the ionized region close to the 
star and the surrounding molecular gas have free-fall den- 
sity profiles, n oc r . At r — r g , the accretion flow 
changes from molecular to ionized. Such a density profile 
wa s indeed found for t he mo lecular envelope of W3 IRS5 
by Ivan der Tak et akl l|2000|) from submillimeter contin- 
uum and line maps. 

The expected flux density of a gravitationally bound 
H II region only depends on the density uq at the ra- 
dius r^where the molecular gas reaches its sound veloc- 
ity (|Ketoll2003^ . Taking T=30 K for the molecular gas, 
r m « 0.35 pc. For d=1.83 kpc, M=20 M and T e =10 4 K, 
the observed flux density of sal mJy at 22-43 GHz is re- 
produced for n wlxlO 5 cm -3 . This estimate agrees to a 
factor of 5 with the valu e of no=2xl0 4 cm -3 derived by 
Ivan der Tak et all l)2000() . We conclude that gravitation 
explains the compactness of the radio sources in W3 IRS5 
which have mid-infrared counterparts. 

4.3. Transient radio sources: proper motions? 

Figure El shows t he positions of radio sour ces A...F de- 
rived b y us and bv lClaussen et al.1 l|l994|) and lWilson et alJ 
(2003J. Sources B and D appear stationary when compar- 
ing the 1989 and 2002 data, but seem to have shifted by 
ss0'/2 in the 1996 data. This shift may be a systematic 
phase error in the 1996 data. The likely cause is an at- 
mospheric 'wedge', or ^100 km-sized parcel of dense air 
partially covering the interferometer and causing a gra- 
dient in atmospheric opacity (M. Reid, priv. comm.). By 
aligning the positions of sources B and D in 1996 with 
those of 1989 and 2002, we derive a shift of Aa = 55±7 
and AS — 146±36 mas. However, since the magnitude of 
the shift is uncertain, we derive proper motion values both 
before and after applying the shift (Table 0}. 

Source C was only detected at two epochs, which is 
sufficient to estimate the magnitude and the direction of 
its motion. For sources A and F, three epochs are avail- 
able, which allows us to solve additionally for the position 
at the first epoch, using the least squares technique. For 
source A, good fits (x 2 /dof ~1) to the a and S motions 
are obtained if no shift is applied. However, for source F, 
the fits are poor (x 2 /dof >10) whether the shift is applied 
or not. One possible explanation for these poor fits are 
deviations from the assumed uniform motions. The 1996 
data thu s may confirm the pr oper motion of component F 
found by I Wilson et al.l l)2003|) and show that components 
A and C may move as well, but do not allow precise mea- 
surement of the magnitude and direction of the motions. 
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Fig. 6. Relative positions of radio sources in W3 IRS5 
before (top) a nd after (bottom) aligning; sources B and D . 
Data are from lClaussen et alJ l|l994h , IWilson et ail l|2003h , 
and Tabled Arrows denote derived proper motions. The 
symbol sizes represent the formal position error of 0'.'06. 



If radio sources A, C and F are internally ionized, 
their ionizing stars must be moving along, since at an 
electron density of 3xl0 6 cm" 3 , the recombination time 
scale is ~1 month while the sources are seen over sev- 
eral years. The free-fall speed in the gravitational poten- 
tial of the molecular g as (M=260 M , i?=60,000 AU: 
Ivan der Tak et aLll2000|) is «3 km s _1 . The potential of 
the star cluster can be estimated by M=20 M Q and 
f?=1000 AU (the typical separation of the radio sources: 
Fig. nj, which gives a free-fall speed of «6 km s _1 . The 
derived proper motions of the radio sources are much 
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Source 


Offset ( 


mas) a 


Proper motion (mas/yr) 




a 


5 


a 


5 






Before Shift 






A 


-0.093 ± 11.550 


0.158 ± 11.223 


5.015 ± 1.801 


8.202 ± 1.624 


C 






-7.813 ± 2.098 


-8.642 ± 2.112 


F 


-2.456 ± 8.277 


-27.419 ± 7.630 


-5.359 ± 1.300 


16.480 ± 0.806 






After Shift 






A 


0.485 ± 11.550 


16.559 ± 11.223 


12.015 ± 1.801 


23.176 ± 1.624 


C 






-0.654 ± 2.098 


10.242 ± 2.112 


F 


0.082 ± 8.277 


26.035 ± 7.630 


1.042 ± 1.300 


16.340 ± 0.806 



a Best fit offset from the measured 1989 position; unavailable for C where only two epochs were measured 



faster than these values, and therefore would not represent 
bound motions. Perhaps these stars were ejected from the 
cluster in a close stellar encounter, and are very young 
runaway OB stars, lik e the BN object in Or ion, which is 
moving at 50 km s" 1 l|Plambeck et alJll995|) . 

We conclude that the evidence for proper motions re- 
mains weak, even with three epochs measured. This shows 
graphically in Fig.|B| the 1996 positions do not generally lie 
between those for 1989 and 2002. Quantitatively, it shows 
in the large error margins in Table ^ The next section 
explores alternative explanations for the transient radio 
sources in W3 IRS 5. 



4.4. Transient radio sources: Shocked clumps? 

The transient radio sources of W3 IRS5 may also be 
explained by shocks which occur when the winds from 
the young O-type stars hit clumps in the surrounding 
molecular material. Such a picture of massive star for- 
mat ion has been descri bed by, e.g., iFranco et al ] lll990h 
and lDvson et al. (2002). Obs ervational supp ort has been 
found in the source Cep A l|Hughesl 12001). which is of 
somewhat lower luminosity and distance than W3 IRS5. 

To estimate the radio emission from wind-shocke d 
clumps, we use the model bv lHollenbach fc McKeel £l989). 
In the limit that the stellar wind is much less dense than 
the molecular clump, the flux Fi of ionizing photons is 
given by 

Fi = n v s F(v s ) 

where no is the density of the clump, v$ the shock ve- 
locity, and F(vs) the fractional ionization of the shocked 
gas. To have F(vs) ~ 1, the shock must fully dissoci- 
ate the H2 clump and heat it to ~10 5 K, so that it 
emits ionizing photons. The required shock velocity is 
~100 km s _1 ; at lower velocities, F(vs) drops exponen- 
tially due to the Boltzmann distribution. Velocities of 
~100 km s _1 are commonly observed for the winds of 
deeply embedded high-mass stars inc luding W3 IRS5 , 
both in hydrogen recombination lines l|Bunn et alj fl995) 



and in CO mid-infrared abso rption lines (jMitchell et alJ 
Il99ltlvan der Tak et allll999|) . 

The emission measure of the shock-ionized clump is 



nil 
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V 10 7 cr 



vs 



m- 3 J V 100km s" 



cm pc 



where 0.3 is the Case B recombination coeffi- 
cient (§ IP|) . The observed values of Z < 100 AU and 
no=10 6 cm -3 agree with this prediction within order of 
magnitude. We conclude that winds from young O-type 
stars shocking and ionizing clumps in the ambient cloud 
provide a viable model for the transient radio emission 
observed in W3 IRS5. 



4.5. Diffuse emission: Envelope structure 

Flux densities for the diffuse mid-infrared emission can 
be obtained by subtracting the point source contributions 
(columns 4, 6 and 10 of TableEJ) from the total flux density 
(column 14). The emission is roughly elliptical in shape, 
with the major axis more or less aligned with the line 
connecting MIR1 and MIR2. Going from short to long 
wavelengths, the axis ratio (measured at the 1% level) 
decreases from wl.6 to wl.0, while the radius (the average 
of the semi-major and semi-minor axes) increases from 
«2300 to «4100 AU. 

The brightness distribution of the diffuse emission is 
consistent with heating by sources MIR1 & MIR2. Short 
wavelengths probe warm dust close to the individual stars, 
so that the emission has two peaks. Longer wavelengths 
probe cooler dust that is far enough away that the dis- 
tances to both stars are about the same, leading to round 
contour shapes. 

At a size of > 1000 AU, the diffuse emission cannot 
be optically thick: even for temperatures as unrealistically 
low as 100 K, the far-infrared luminosity is exceeded, even 
with zero foreground extinction. The envelope must also 
have a low mid-infrared optical depth to give us a view of 
the central objects. 

The envelope o f W3 IRS5 was modeled by 
ICamnbell et al based on far-infrared data, and by 
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Ivan der Tak et all ij^OOO). based on submillimeter data. 
Assuming single power laws for the density structure, 
these models have r(100 /mi)~l, or r(10 ^m)^40, much 
higher than estimated above. The W3 IRS5 core is em- 
bedded in a large-scale molecular cloud, and the observed 
submillimeter emission could contain contibutions from 
the background cloud, but not more than 50%. Neither 
the mid-infrared nor the submillimeter images suggest de- 
viations from spherical symmetry stronger than a modest 
flattening (axis ratio <2). 

To reconcile the submillimeter and far-infrared data 
with the mid-infrared data, one may explore broken power 
laws for the density distribution, or a combination of dense 
shells and power laws. 

5. Conclusions and Outlook 

Observations at subarcsecond resolution at mid-infrared 
and radio wavelengths have led to a detailed picture of 
the W3 IRS5 region, shown schematically in FigureEJand 
described below. 

— The two bright mid-infrared sources with radio emis- 
sion probably are deeply embedded high-mass stars. 
They are both close to groups of H 2 masers. The 
measured mid-infrared diameters are consistent with 
blackbody emission at T=390 K and providing all of 
the far-infrared luminosity of W3 IRS5. 

— A third, weaker mid-infrared and radio source with 
associated H2O masers is probably a somewhat later- 
type star which is energetically unimportant. 

— Radio source A has an optically thick radio spectrum, 
and may have a counterpart at long mid-infrared wave- 
lengths (J 17.65 /mi). It may be an extremely deeply 
embedded high-mass star. The three power sources 
of W3 IRS5 t hen have L»40,000 e ach, which in 
the models of iMaeder fc Mevnetl 1^.989) makes them 
«20 Mq stars (ZAMS spectral type 08). 

— The region shows several transient radio sources. Some 
of these may represent runaway OB stars, but most 
are probably clumps in the ambient material which 
are ionized and destroyed by shocks with the winds of 
the O-type stars. 

— The low silicate optical depth suggests that no under- 
lying silicate emission is present. This is most easily 
explained by a cavity separating the high-mass stars 
from their envelope. Perhaps the cavity was blown by 
the slow spherical outflow traced by the H2O masers. 

— The far-infrared and submillimeter emission, as well 
as the low-velocity CO mid-infrared absorption, arise 
in the large-scale envelope. The dense stellar cluster 
visible in the near-infrared (and X-ray) is embedded 
in the same envelope. 

In the future, subarcsecond monitoring of W3 IRS5 
at high radio frequencies (VLA-A) is necessary to test the 
'proper motion' and 'shocked clump' hypotheses. If proper 
motions are confirmed, the shapes of the orbits will be a 



test of the 'runaway star' hypothesis, and will constrain 
the dynamics of this young cluster. The emission mech- 
anism should be studied by simultaneous observations at 
three or more wavelengths. If the e-VLA does not provide 
the sensitivity necessary to do this, ALMA will. 

Mid-infrared imaging at wavelengths > 20 /im at sub- 
arcsecond resolution is necessary to search for a mid- 
infrared counterpart to radio source Q2=K3=A. The 
main requirements are higher sensitivity and dynamic 
range than was achieved here; the higher angular reso- 
lution offered by MIDI on the VLTI will be more use- 
ful to search for fine structure. Spatially resolved mid- 
infrared spectroscopy is necessa ry to assign the h igh- 
velocity CO absorption features IjMitchell et alJll991f) to 
particular stellar components. This may be a good project 
for VLT/CRIRES. Future radiative transfer modeling ef- 
forts should consider broken power laws for the density 
distribution in the envelope of W3 IRS 5. 
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